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ABSTRACT 

The lensing magnification effect due to large-scale structure is statistically measurable by cor- 
relation of size fluctuations in distant galaxy images as well as by cross-correlation between 
foreground galaxies and background sources such as the QSO-galaxy cross-correlation. We 
use the halo model formulation of Takada & Jain (2003) to compute these magnification- 
induced correlations without employing the weak lensing approximation, \i w 1 + 2k. Our 
predictions thus include the full contribution from non-linear magnification, 6fi £ 1, that is 
due to lensing halos. We compare the model prediction with ray-tracing simulations and find 
excellent agreement over a range of angular scales we consider (0'.5 < 9 < 30')- In addition, 
we derive the dependence of the correlation amplitude on the maximum magnification cutoff 
Mmax, which is necessary to introduce in order to avoid the contributions from strong lens- 
ing events. For a general correlation function parameterized as (/i p /) (/ is any cosmic field 
correlated with the magnification field), the amplitude remains finite for p < 1 and diverges 
for p > 1 as /i max — > oo, independent of details of the lensing mass distribution and of the 
separation angle. This consequence is verified by the halo model as well as by the simulations. 
Thus the magnification correlation with p < 1 has a practical advantage in that it is insensitive 
to a selection effect of how strong lensing events with fi 1 are observationally excluded 
from the sample. 

The non-linear magnification contribution enhances the amplitude of the magnification 
correlation relative to the weak lensing approximation, and the non-linear correction is more 
significant on smaller angular scales and for sources at higher redshifts. The enhancement 
amounts to 10 — 25% on arcminite scales for the QSO-galaxy cross-correlation, even after in- 
clusion of a realistic model of galaxy clustering within the host halo. Therefore, it is necessary 
to account for the non-linear contribution in theoretical models in order to make an unbiased, 
cosmological interpretation of the precise measurements expected from forthcoming massive 
surveys. 
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1 INTRODUCTION 

Gravitational lensing caused by the large-scale structure is now rec- 
ognized as a powerful cosmological tool (Mellier 1999; Bartelmann 
& Schneider 2001 for thorough reviews). The gravitational deflec- 
tion of light causes an increase or decrease of the area of a given 
patch on the sky depending on which the light ray passes prefer- 
entially through the overdense or underdense region. Accordingly, 
this causes an observed image of source to be magnified or de- 
magnified relative to the unlensed image, since lensing conserves 
the surface brightness and the received luminosity is proportional to 
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the solid angle of the image. Large magnifications are observed in a 
strong lensing system that accompanies multiple images or largely 
deformed images. It has been also proposed that mild or weak mag- 
nifications is measurable in a statistical sense. The magnification 
leads to an enhancement in the flux-limited number counts of back- 
ground sources around foreground sample that traces the lensing 
mass distribution. Based on this idea, numerous works have inves- 
tigated the QSO-galaxy cross correlation theoretically (e.g., Broad- 
hurst, Taylor & Peacock 1995; Bartelmann 1995) as well as obser- 
vationally (e.g., Benitez & Martinez-Gonzalez 1997; Benitez et al. 
2001; Gaztanaga 2003; also see Bartelmann & Schneider 2001 for 
a thorough review). Further, Jain (2002) recently proposed that the 
magnification effect can be extracted by statistically dealing with 
size fluctuations of distant galaxy images, though it is a great chal- 
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lenging for existing data yet. Forthcoming massive surveys such 
as SDSS 1 , DLS 2 , the CFHT Legacy survey 3 as well as future pre- 
cise imaging surveys such as SNAP 4 , Pan-STARRS 5 and LSST 6 
allow measurements of these magnification effects at high signifi- 
cance. Therefore, it is of great importance to explore how this type 
of method can be a useful cosmological tool and complementary to 
the established cosmic shear method, which measures correlation 
between lensing induced ellipticities of distant galaxy images (e.g., 
Hamana et al. 2002 and Jarvis et al. 2003 and references therein). 

The magnification field fi(8) in a given direction on the sky 
is expressed (e.g, Schneider, Falco & Ehlers 1992) as 

^)=i(i-^)) 2 -7 2 wr 1 . CD 

Here k and 7 are the convergence and shear fields, which are fully 
determined by the mass distribution along the line of sight. This 
equation shows the nonlinear relation between the magnification 
and the convergence and shear, and indicates that the magnification 
increases with k and 7 very rapidly and becomes even formally 
infinite when the lensing fields k, 7 ~ O(l). However, as long 
as we are concerned with the magnification related statistics due 
to the large-scale structure, strong lensing events (/1 ^> 1) should 
be removed from the sample to prevent the large statistical scat- 
ters. This will be straightforward to implement, if the strong lens- 
ing accompanies multiple images or largely deformed images. On 
the other hand, modest magnification events (<5/i £ 1) make it rel- 
atively difficult to identify and are likely included in the sample for 
the blind analysis, since the magnification is not a direct observable. 
Therefore, the magnification statistics rather requires a more care- 
ful study of the selection effect than the cosmic shear (e.g., Barber 
& Taylor 2003), which we will carefully address. 

The simplest statistical quantity most widely used in cosmol- 
ogy is the two-point correlation function (2PCF). For our purpose, 
the magnification field is taken as either or both of the two fields 
entering into the correlation function. However, it is not straight- 
forward to analytically compute the magnification 2PCF because of 
the non-linear relation between /i and n, 7, where the latter fields 
are easier to compute in a statistical sense based on a model of 
the mass power spectrum. For this reason, the conventional method 
of the magnification statistics employs the weak lensing approx- 
imation /i ~ 1 + 2k (e.g., Bartelmann 1995; Dolag & Bartel- 
mann 1997; Sanz et al. 1997; Benitez & Martinez-Gonzalez 1997; 
Moessner & Jain 1998; Benitez et al. 2001; Menard & Bartelmann 
2002; Gaztanaga 2003; Jain et al. 2003). However, it is obvious that 
this type of method is valid only in the limit k, 7 <JC 1 and likely 
degrades the model accuracy on non-linear small scales. In fact, 
using the ray-tracing simulations Menard et al. (2003) clarified the 
importance of the non-linear magnification contribution to the mag- 
nification statistics (see also Barber & Taylor 2003). It was shown 
that the perturbative treatment breaks down over a range of angu- 
lar scales of our interest. Although a promising method to resolve 
this issue is to employ ray-tracing simulations, to perform multiple 
evaluations in model parameter space requires sufficient number of 
simulation runs, which is relatively prohibitive. 

Therefore, the main purpose of this paper is to develop an 
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analytic method to compute the magnification induced correlation 
function without employing the weak lensing approximation. To do 
this, we use the halo model to describe gravitational clustering in 
the large-scale structure, following the method developed in Takada 
& Jain (2003a,b,c, hereafter TJ03a,b,c). The model prediction of 
the QSO-galaxy cross-correlation is also developed by incorporat- 
ing a realistic model of galaxy clustering within the host halo into 
the halo model. Although the halo model rather relies on the sim- 
plified assumptions, the encouraging results revealed so far are that 
it has led to consistent predictions to interpret observational results 
of galaxy clustering as well as to reproduce simulation results (e.g., 
see Seljak 2000; Zehavi et al. 2003; Takada & Jain 2002, hereafter 
TJ02; TJ03a,b,c; also see Cooray & Sheth 2002 for a review). 

Another purpose of this paper is to explore how the magnifi- 
cation related statistics can probe the halo structure. The non-linear 
magnifications arise when the light ray emitted from a source en- 
counters an intervening mass concentration, i.e., dark matter halo 
such as galaxy or cluster of galaxies. It is known that strong lensing 
of \x 3> 1 can be used to probe detailed mass distribution within a 
halo (e.g., Hattori et al. 1999). Similarly, modest non-linear magni- 
fications of <5/i ^ 1 could lead to a sensitivity of the magnification 
statistics to the halo structure in a statistical sense, as investigated in 
this paper. A fundamental result of cold dark matter (CDM) model 
simulations is that the density profiles of halos are universal across 
a wide range of mass scales (e.g., Navarro, Frenk & White 1997, 
hereafter NFW). On the other hand, some alternative models such 
as self-interacting dark matter scenario (Spergel & Steinhardt 2000) 
have been proposed in order to reconcile the possible conflicts be- 
tween the simulation prediction and the observation. If dark matter 
particle has a non-negligible self-interaction between themselves, 
the effect is likely to yield a drastic change on the halo profile 
compared to the CDM model prediction (Yoshida et al. 2000). The 
halo structure thus reflects the dark matter nature as well as detailed 
history of non-linear gravitational clustering. Hence, exploring the 
halo profile properties with gravitational lensing can be a direct test 
of the CDM paradigm on scales £ Mpc, which is not attainable 
with the cosmic microwave background measurement. 

The plan of this paper is as follows. §2 presents the basic 
equations relevant for cosmological gravitational lensing and then 
briefly summarize two promising methods to statistically measure 
the magnification effect. In §3 we develop an analytic method to 
compute the magnification correlation functions based on the halo 
model. In §4, we derive an asymptotic behavior of the correla- 
tion amplitude for large magnifications. In §5 we qualitatively test 
the halo model prediction and the asymptotic behavior using ray- 
tracing simulations. The realistic model of the QSO-galaxy cross- 
correlations is also derived. §6 is devoted to a summary and dis- 
cussion. Throughout this paper, without being explicitly stated, we 
consider the ACDM model with Q. m o = 0.3, Q.\o = 0.7, flbu, 
h = 0.7 and as = 0.9. Here ftmo, Qbo and fl\o are the present- 
day density parameters of matter, baryons and the cosmological 
constant, h is the Hubble parameter, and as is the rms mass fluctu- 
ation in a sphere of radius 8/i _1 Mpc. 



2 PRELIMINARIES 

2.1 Magnification of gravitational lens 

The gravitational deflection of light ray induces a mapping between 
the two-dimensional source plane (S) and the image plane (I) (e.g., 
Schneider, Ehlers & Falco 1992): 
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<Bf = A,* 1 /, (2) 

where <5#i is the separation vector between points on the respec- 
tive planes. The lensing distortion of an image is described by the 
Jacobian matrix Aij defined as 

A ~ ( 1 " K " 71 ~ 72 ^ (3) 

where k is the lensing convergence and 71 and 72 denote the tidal 
shear fields, which correspond to elongation or compression along 
or at 45° to x-axis, respectively, in the given Cartesian coordinate 
on the sky. The k and ji depend on angular position, although we 
have omitted showing it in the argument for simplicity. Since the 
gravitational lensing conserves surface brightness from a source, 
the lensing magnification, the ratio of the flux observed from the 
image to that from the unlensed source, is given by determination 
of the deformation matrix, yielding equation (1). In the weak lens- 
ing limit k, 7 <§C 1, one can Taylor expand the magnification field 
as /j, ss 1 + 2k, as conventionally employed in the literature to com- 
pute the magnification statistics. The weak lensing approximation 
ceases to be accurate as k, 7 — > O(l). For example, k = 0.5 (and 
simply 7 = 0) leads to factor 2 difference as /j = 4 and 1 + 2k = 2. 

In the context of cosmological gravitational lensing, the con- 
vergence field is expressed as a weighted projection of the three- 
dimensional density fluctuation field between source and observer: 

k(0)= r H d X W(x)S[x,d A (x)0], (4) 
Jo 

where x ls the comoving distance, (Ia(x) is m e comoving angu- 
lar diameter distance to x> an d X h is the distance to the Hub- 
ble horizon. Note that X is related to redshift 2 via the relation 
dx = dz/H(z) (H(z) is the Hubble parameter at epoch z). 
Following the pioneering work done by Blandford et al. (1991), 
Miralda-Escude (1991) and Kaiser (1992), we used the Born ap- 
proximation, where the convergence field is computed along the 
unperturbed path, neglecting higher order terms that arise from cou- 
pling between two or more lenses at different redshifts. Using ray- 
tracing simulations of the lensing fields, Jain et al. (2000) proved 
that the Born approximation is a good approximation for lensing 
statistics (see also Van Waerbeke et al. 2001; Vale & White 2003). 
The lensing weight function W is given by 

W(x) = ^ m0 Hla-\ X )d A { X ) 

where f s ( X s) is the redshift selection function normalized as 
j* H dxfs{x) = 1- In this paper we assume all sources are at 
a single redshift z s for simplicity; f s ( X ) = $r>(x — X«)- ^0 is 
the Hubble constant (H = lOO/i km s _1 Mpc _1 ). Similarly, the 
shear fields are derivable from the density fluctuation fields, but 
the relation is non-local due to the nature of the gravitational tidal 
force. In Fourier space, we have the simple relation between k and 
7: 

7i(Z) = cos2ipik(l), 72(0 = sin 2</?jre(Z), (6) 
where the Fourier-mode vector is I = /(cos sin tp}). 



2.2 Methodology for measuring the magnification statistics 

There are two promising ways for measuring the lensing magnifi- 
cation effect statistically, which are likely feasible for forthcoming 
and future surveys. Here we briefly summarize the methodology. 

Gravitational magnification has two effects. First, it causes the 
area of a given patch on the sky to increase, thus tending to dilute 
the number density observed. Second, sources fainter than the lim- 
iting magnitude are brightened and may be included in the sample. 
The net effect, known as the magnification bias, depends on how 
the loss of sources due to dilution is balanced by the gain of sources 
due to flux magnification. Therefore, it depends on the slope s of 
the unlensed number counts of sources No(m) in a sample with 
limiting magnitude m, s = d log No(m) / dm. Magnification by 
amount /j changes the number counts to (e.g., Broadhurst, Taylor 
& Peacock 1995; Bartelmann 1995); 

N'(m) = N {m)n 2 - 5s - 1 . (7) 

For the critical value s = 0.4, magnification does not change the 
number density; it leads to an excess for s > 0.4, and a deficit 
for s < 0.4. Let n\{0) be the number density of foreground 
sample with mean redshift (2)1, observed in the direction 6 on 
the sky, and 712(6) that of the source sample with a higher mean 
redshift (2)2 > {z)i. Thus, even if there is no intrinsic correla- 
tion between the two populations, magnification induces the non- 
vanishing cross-correlation: 

m = («5m(0i)<5n 2 (0 2 )) |0i _ 02|=e 

= (^iMM^)] 2 - 5 *- 1 ),^.^, (8) 

where <5n;(0) = (rii(G) — fii))/ni, with fn the average number 
density of the ith sample. Here (• • ■) denotes the ensemble average 
and observationally means the average over all pairs separated by 8 
on the sky. Based on this idea, numerous studies have confirmed the 
existence of the enhancement of the QSO number counts around 
foreground galaxies, i.e., the QSO-galaxy cross correlation (e.g., 
Benitez et al. 2001, Gaztanaga 2003 and references therein; also 
see Bartelmann & Schneider 2001 for a review). Although these 
results are in qualitatively agreement with the magnification bias, 
in most cases the amplitude of the correlation is much higher than 
that expected from gravitational lensing models. The excess might 
be due to the non-linear magnification contribution from massive 
halos (e.g., see the discussion around Figure 28 in Bartelmann & 
Schneider 2001). However, to make such a statement with confi- 
dence, it is necessary to further explore an obstacle in the theoret- 
ical model, the bias relation between the galaxy and mass distri- 
butions. This still remains uncertain observationally and theoreti- 
cally. In particular, on small scales ( ^ Mpc), it is crucial to model 
how to populate halos of given mass with galaxies, known as the 
halo occupation number (e.g., Seljak 2000; TJ03b). Recently, Jain, 
Scranton & Sheth (2003) carefully examined the effect of the halo 
occupation number on the magnification bias and showed that the 
parameters used in it yield a strong sensitivity to the predicted cor- 
relation. For example, possible modification for types of galaxies 
leads to a change of a factor of 2-10 in the expected signal on ar- 
cminute scales. Further, we will show that the non-linear magnifi- 
cation correction is also important to make the accurate prediction. 

Second, Jain (2002) recently proposed a new method for mea- 
suring the magnification effect of the large-scale structure, based 
on the work of Bartelmann & Narayan (1995). The lensing effect is 
to increase or decrease an observed image of galaxy relative to the 
unlensed image, depending on which the light ray travels prefer- 
entially through overdense or underdense region that corresponds 
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to n > 1 or n < 1. That is, the area 5* and characteristic radius 
-R(oc S 1 / 2 ) are changed by the magnification field p as 



S — > S/i, 



1/2 



(9) 



Although the unlensed image is not observable, this effect can be 
statistically extracted as follows. We can first estimate the mean 
size of source galaxies from the average over the sample avail- 
able from a given survey area, under the assumption (pi) = 1 or 
(pi 1 / 2 ) = 1. Then, the two-point correlation function (2PCF) of the 
size fluctuations can be computed from the average over all pairs 
of galaxies separated by the angle considered, in analogy with the 
2PCF of the cosmic shear fields (the variance method was consid- 
ered in Jain 2002). The reason that this method is not yet feasible 
is that it requires a well controlled estimate of the unlensed size 
distribution as well as systematics (photometric calibration, reso- 
lution for sizes and PSF anisotropy). Hence, space based imaging 
surveys will make possible the measurement of galaxy sizes with 
a sufficient accuracy hard to achieve from the ground so far. This 
method can potentially be a precise cosmological probe as the cos- 
mic shear measurement, because it is free from the galaxy bias un- 
certainty. Further, the non-linear nature of the magnification could 
lead to complementarity to the cosmic shear measurement, as we 
will discuss below. 



3 HALO APPROACH TO MAGNIFICATION STATISTICS 
3.1 Halo profile and mass function 

We use the the halo model of gravitational clustering to compute 
the magnification statistics, following the method in TJ03a,b,c. Key 
model ingredients are the halo profile ph{r), the mass function 
n(M), and the halo bias b(M), each of which is well investigated 
in the literature. 

As for the halo profile, we employ an NFW profile truncated 
at radius riso, which is defined so that the mean density enclosed 
by sphere with riso is 180 times the background density. Within a 
framework of the halo model, we need to express the NFW profile 
in terms of Miso and redshift z. To do this, we first express the 
two parameters of the NFW profile, the central density parameter 
and the scale radius, in terms of the virial mass and redshift, based 
on the spherical top-hat collapse model and the halo concentration 
parameter of Bullock et al. (2001) (see TJ03b,c for more details). 
Then, following Hu & Kravtsov (2002), we employ a conversion 
method between the virial mass and Miso in order to re-express 
the NFW profile in terms of Miso and z. In what follows, we will 
often refer halo massed M as Miso for simplicity. 

For the mass function, we employ the Sheth-Tormen mass 
function (Sheth & Tormen 1999): 



n(M)dM = 



M 
Po 
M 



f(y)dv 
A[l + [av) 



auexp 



I av\ dv 



where v is the peak height defined by 

d. ( : ) 

v - 



D(z)a{M) 



(10) 



(11) 



Here <r(M) is the present-day rms fluctuation in the mass density, 
smoothed with a top-hat filter of radius Rm = (3M/47rpo) 1 ^ 3 , S c 
is the threshold overdensity for the spherical collapse model and 
D(z) is the linear growth factor (e.g., Peebles 1980). The numer- 
ical coefficients a and p are taken from the results of Table 2 in 



White (2002) as a = 0.67 and p = 0.3, which are different from 
the original values a = 0.707 and p = 0.3 in Sheth & Tormen 
(1999). Note that the normalization coefficient A = 0.129. The 
main reason we employ the halo boundary riso is that the mass 
function measured in iV-body simulations can be better fitted by the 
the universal form (10) when one employs the halo mass estimator 
of Miso, than using the virial mass estimator, as carefully exam- 
ined in White (2002). To maintain the consistency, we also employ 
the halo biasing b(M) given in Sheth & Tormen (1999) with the 
same a and p parameters, which is needed for the 2-halo term cal- 
culation. Note that J dM b(M)n(M)M/p ~ 1 for the halo bias 
model and the mass function we employ (e.g., Seljak 2000). 

3.2 Lensing fields for an NFW profile 

For an NFW profile, we can derive analytic expressions to give the 
radial profiles of convergence, shear and magnification. 

For a given source at redshift z s , the convergence profile for a 
halo of mass M at z, denoted by km, can be given as a function of 
the projected radius from the halo center: 



<Ia(Xs 



(12) 



where r s is the scale radius, 8 S its projected angular scale and / = 
l/[ln(l + C180) — Ci8o/(l + C180)] (see below for the definition 
of ci8o). The explicit form of F(x) is given by equation (27) in 
TJ03b. For an axially symmetric profile, the shear amplitude can 
be derived as 



7m(6>) 



R(6) - k{6) 

- 1 d A (x)dA(Xs~x) Mf 



AnGa 



<Ia{Xs 



2nr 



(13) 



where G(x) is given by equation (16) in TJ03c. These expressions 
of km and j M differ from those given in Bartelmann (1996), which 
are derived from an infinite line-of-size projection of the NFW pro- 
file under the implicit assumption that the profile is valid (even far) 
beyond the virial radius. TJ03b,c carefully verified that employ- 
ing the expressions (12) and (13) in the real-space halo model is 
essential to achieve the model accuracy as well as the consistency 
with the Fourier-space halo model well investigated in the literature 
(e.g., Seljak 2000). 

Throughout this paper we employ the halo boundary riso, not 
the virial radius, as stated in §3.1. For this setting, we have to re- 
place parameter c used in F(x) and G(x) in TJ03b,c with the ratio 
of the scale radius to riso, ciso- The relation between ciso and c 
is given by ciso = criso /r v i r - Note that, even if we use the virial 
boundary, the results shown in this paper are almost unchanged. 

Given the convergence and shear profiles for a halo of mass 
M, the magnification profile is given by [xm (0) = |(1— km{0)) 2 — 
7m This equation implies that [im becomes formally infin- 

ity at some critical radius when the denominator becomes zero. The 
radius in the lens plane is called the critical curve, while the corre- 
sponding curve in the source plane is the caustic curve. 

Any lensing NFW halo inevitably provides finite critical ra- 
dius, if we have an ideal angular resolution. This is because the 
convergence km varies from zero to infinity with changing radius 
from 6<i8o to 0, while the shear remains finite over the range (see 
Figure 2 in TJ03c). Figure 1 plots the tangential critical radius for 
a lensing NFW halo against the mass. Note that an NFW profile 
produces two critical radii - an outer curve causes a tangentially 
deformed image around it, while the inner one causes a radially de- 
formed image. The solid and dashed curves are the results for halos 
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(19) and (20) in TJ03c): 



M [M e ] 



Figure 1. The tangential critical radius of NFW profile against halo masses. 
The solid and dashed curves are the results for lensing halos at redshifts of 
zi = 0.4 and 0.8 for source redshifts of z s = 1 and z s = 3, respectively. 
Even if we consider a lower lens redshift, the result does not largely change, 
as shown by the dot-dashed curve. 



at redshifts zi — 0.4 and 0.8 for source redshifts of z s — 1.0 and 
3.0, respectively. One can see that the critical radius has a strong 
dependence on halo masses and is larger for z s — 3 than for z s = 1 
due to the greater lensing efficiency. These critical curves do not 
largely change even if we consider a lower lens redshift than the 
peak redshift, as shown by the dot-dashed curve. Even massive ha- 
los with M ~ 1O 15 M provide the critical radii of < O'.l. The 
scale is below relevant angular scales for the magnification statis- 
tics of our interest. However, this does not mean that the non-linear 
magnification correction to the correlation function appears only 
on scales ,$ 1'. Rather, modest non-linear magnifications {p, £ 2) 
lead to the strong impact, since such magnifications have larger 
cross section, as will be shown in detail. 



3.3 Real-space halo model approach 

In the following, we construct the halo model method to compute 
the magnification induced correlation function. First, we simply 
consider the 2PCF defined as 



(,(«) = {^)-iM0 + fl)-i]), 



(14) 



where fj, — 1 is the magnification fluctuation field. This 2PCF is 
observable from size fluctuations of distant galaxy images, as dis- 
cussed in §2.2. 

From a picture of the halo model, can be expressed as 

a sum of correlations between the magnifications fields within a 
single halo (1-halo term) and between two different halos (2-halo 
term): 



(15) 



It is straightforward to extend the real-space halo model developed 
in TJ03a,b,c to compute the 1-halo term contribution, which has 
dominant contribution on small non-linear scales (see equations 



XH ,2y 



x d<j> #0[mm(c/>)-1][pm(|0 + 0|)-1],(16) 
Jo Jo 

where we have introduced the polar coordinate d 2 cf> — <f>d(j)dip, 
d 2 V/d X dn = d\( X ) for a flat universe and we can set the sepa- 
ration vector 6 to be along the first axis from statistical symmetry, 
thereby \<p + 0\ = (4> 2 + 6 2 + 2(f>6 cos f) 1/2 . We have assumed a 
uniform distribution of sources on the sky and ignored a probabil- 
ity of multiple images and an increase or decrease in sampling of 
the images due to the lensing itself (which is a higher order correc- 
tion and can be safely neglected as shown in Hamana 2001). The 
equation above implies that the 1-halo term contribution is given 
by summing lensing contributions due to halos along the line-of- 
sight weighted with the halo number density on the light cone. Note 
that the integration range of d 2 <p is taken as the infinite area, tak- 
ing into account the non-local property of the shear field that is 
non- vanishing at radius outside the halo boundary. In practice, set- 
ting the upper bound of Jd(j> to be three times the projected radius 
#180 gives the same result, to within a few percents. The validity of 
the real-space halo model formulation was carefully investigated in 
TJ03b,c. 

As discussed in §3.2 and shown in Figure 1, an NFW pro- 
file always provides finite critical curves, where the magnification 
formally becomes infinity. Therefore, to make the halo model pre- 
diction, we introduce a magnification cutoff /i max in the calculation 
- the integration range of fd 2 (fi is confined to the region satisfying 
the condition [i u < (i ma , for a given halo. This procedure is some- 
how similar to what is done in the measurement from ray-tracing 
simulations, where a masking of high magnification events is em- 
ployed to avoid a significant statistical scatter (e.g., see Menard et 
al. 2003; Barber & Taylor 2002). Thus the halo model allows for a 
fair comparison of the prediction with the simulation result. How- 
ever, note that the procedure taken ignores the lensing projection 
effect for simplicity: exactly speaking, the magnification should be 
given by the lensing fields between source and observer, not by 
those of individual halo. 

Similarly, based on the real-space halo model, the 2-halo term 
can be expressed as 



l3 x d x dn 



dMn(M) dM' n(M') 



xb(M)b(M') 



x d 2 4> \n M (<f>) - 1] d 2 cj>' [uttW) - 1] 



x) J {i\e-<j>-<t>' 



2iy \ d A 



(17) 



where Jo(x) is the zero-th order Bessel function and P L (k) is the 
linear mass power spectrum at epoch \ as given by P L (k; \) = 
D 2 (x)P(k; t ). The term including P L (k) in the third line on the 
r.h.s denotes the angular two-point correlation function between 
different halos of masses M and M', which is derived using Lim- 
ber's equation and the flat-sky approximation (e.g., Blandford et al. 
1991; Miralda-Escude 1991; Kaiser 1992). We similarly introduce 
the maximum magnification cutoff in the 2-halo term calculation 
as in the 1-halo term. The equation above means that we have to 
perform an 8-dimensional integration to get the 2-halo term, which 
is computationally intractable. For this reason, we employ an ap- 
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proximated way to be valid when angular separation between the 
two halos taken is much greater than their angular virial radii. This 
leads to a simplified expression of the 2-halo term: 



C h (0) 



I 

Jo 



dx 



d 2 V 



dxdQ, 

JdMb{M)n{M) J<t>d<j) 2n(n M (4>)-l) 

x) Mm- 



2% 



k - 



dAixY 



(18) 



This allows us to get the 2-halo term by a three-dimensional inte- 
gration, because the integrations of the second and third lines on the 
r.h.s can be done separately before performing the x-integration. It 
is worth checking the consistency of the 2-halo term above with 
the limiting case £ M w 4£ K , which is valid in the weak lensing limit 
/i Ri 1 + 2k for k, 7 <C 1. In this limit, the term contained in the 
second line on the r.h.s of equation (18) can be rewritten as 

JdM b(M)n{M) Jd 2 (b(nM{<l>)-l) 



dMb(M)n(M) d 2 <t>2K M ((f>) 



2W( X )d 



a\x) j< 



M 

dM b{M)n{M) — , 
Po 



(19) 



where the second equality is derived from equations (25) and (28) 
in TJ03b. Therefore, substituting this result into the 2-halo term 
(18) yields 



4 / d X W\ X )d-/{x) 
M 



dM b(M)n(M) 



f 



Po 

dkr x ) MW) 



(20) 



The 2-halo term (18) is thus reduced to four times the 2-halo term 
of the convergence 2PCF in the weak lensing limit, as expected. 
However, the important point of the 2-halo term (18) is that it can 
correctly account for the contribution of non-linear magnifications 
(/i £ 2) on large scales. 

Replacing P 2h (k) in equation (20) with a model of the non- 
linear mass power spectrum leads to the conventional method for 
predicting 4£ K (0) as an estimator of as employed in the lit- 

erature (e.g., Dolag & Bartelmann 1997; Sanz et al. 1997; Jain et 
al. 2003). 



4 ASYMPTOTIC BEHAVIOR OF MAGNIFICATION 
STATISTICS FOR HIGH MAGNIFICATIONS 

As stated in the preceding section, we need to introduce a maxi- 
mum magnification cutoff in the model prediction to avoid the con- 
tribution from a formally emerged infinite magnification. In prac- 
tice, strong lensing events of /i 3> 1, which are identified by mul- 
tiple images or largely deformed image, can be removed from the 
sample of the magnification statistics. However, without the clear 
signatures, it is hard to make the distinct selection and therefore 
modest magnification events (5fi £ 1) are likely included in the 



sample for the blind analysis, because the magnification is not a di- 
rect observable. In the following, we clarify how the magnification 
statistics depend upon large magnification events (/i ^> 1). 

Meanwhile, we restrict our discussion to point-like sources 
for simplicity. High magnifications (fj, >• 1) arise from images 
in the vicinity of the critical curve that is caused by an intervening 
mass concentration, such as halos. For any finite lens mass distri- 
bution, the critical curve must form a closed non-self-intersecting 
loop. Based on the catastrophe theory, it was shown in Blandford 
& Narayan (1986; also see Chapter 6 in Schneider et al. 1992) that 
the magnification of an image at a perpendicular distance AO from 
the (fold) critical curve scales asymptotically as 

1 



jX oc 



AO' 



(21) 



This argument holds independently of details of the lensing mass 
distribution, although the proportional coefficient does depend on 
the mass distribution. 

To keep generality of our discussion, let us consider a corre- 
lation function between the magnification field and some cosmic 
fields expressed as 

mP) = (HOi)] P f(02)) l0l _ 02l=e , (22) 

where p is an arbitrary number. The field / is allowed to be consti- 
tuted from any cosmic fields which are correlated with the magnifi- 
cation field. Therefore, the following argument holds for high-order 
moments beyond the two-point correlations if one takes products of 
the cosmic fields for /, e.g. / = 8(62)8(63). 

Suppose that an intervening halo provides the critical curve in 
the lens plane for a given source redshift, as this is the case for an 
NFW halo (see Figure 1). Then, let us consider how high magnifi- 
cations in the vicinity of this critical curve contribute to the mag- 
nification statistics. By introducing an upper bound on the magni- 
fication fields as /i < /i max , we can address how the magnification 
statistics depend on the cutoff /i max and what is the asymptotic be- 
havior for the limiting case /i max — > 00. From equation (21), the 
cutoff /imax corresponds to a lower limit on the perpendicular dis- 
tance from the critical curve, say AO > e (e — > corresponds to 
/imax — » 00). As can be seen from equations (16) and (18), a pic- 
ture of the halo model leads us to compute the high magnification 
contribution to the magnification-induced correlation like equation 
(22) by the integration 

£(0;p,/w)~ f d 2 ST ^f(\s + e\), (23) 



|As|p" 



where As is the perpendicular distance from the critical curve and 
the integration range is confined to the area dS c subject to the con- 
dition /i < /imax- The integration (23) results in one-dimensional 
integration for high magnifications around the critical curve, in 
analogy with equation (5.16) in Blandford & Narayan (1986) to 
derive the asymptotic, integral cross section for the strong lens- 
ing events that produce multiple images 7 . Hence, the leading order 
contribution of /i ma x can be expressed as 



£(0;P,Mmax) 



1-p 

max 1 



ln(/w 

Pmax 1 



P < 1 
p = i 
p > 1. 



(24) 



7 The discussion in this paper as well as in Blandford & Narayan (1986) 
employs the assumption that asymptotic dependence of the magnification 
statistics on high magnifications is mainly due to images around the fold 
caustics and therefore ignores the contribution from the cusp caustic. This 
is likely to be a good approximation as shown in Mao (1992). 
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where we have assumed that variation in the field / does not largely 
change for the relevant integration range. The equation above leads 
to an intriguing consequence: the amplitude of the magnification 
correlation is finite for p < 1 for the limiting case (i ma , — > oo, 
while it diverges for p > 1. Thus, the statistics with p < 1 is 
practically advantageous in that it is insensitive to the uncertainty of 
which magnification cutoff should be imposed for a given sample. 
Furthermore, the asymptotic behavior does not explicitly depend 
on the separation angle 8 and therefore it holds even for large 6. 
This means that the divergence £ — > oo for p > 1 formally occurs 
even on degree scales, which is opposed to a naive expectation that 
the weak lensing approximation is safely valid on these scales. It is 
also worth stressing that this behavior is expected to hold for any 
lensing mass distribution once the critical curve appears, although 
the proportionality coefficient of £ should depend on details of the 
mass distribution. These will be quantitatively tested by the halo 
model prediction as well as by the ray-tracing simulation. 

In reality, a finite source size imposes a maximum cutoff on 
the observed magnification and thus the infinite magnification does 
not happen, even if a source sits on the caustic curve in the source 
plane (see Chapters 6 and 7 in Schneider et al. 1992; Peacock 1982; 
Blandford & Narayan 1996). For example, if the source is a circular 
disk with radius rs and uniform surface brightness, the maximum 
magnification is given by ^ m ax oc 1/ *Jr~s. Therefore, to develop an 
accurate model prediction requires a knowledge of unlensed source 
properties such as the size and the surface brightness distribution in 
addition to modeling the lensing mass distribution. In particular, 
this could be crucial if one considers the magnification statistics 
(22) with p > 1, since it is sensitive to large magnification events. 



5 RESULTS 

5.1 Ray-tracing simulations 

To test the analytic method developed in §3, we employ the ray- 
tracing simulations. We will use the simulation for the current con- 
cordance ACDM model with f2 m0 = 0.3, Q\ = 0.7, fi b o = 0.04, 
h = 0.7 and <T8 = 0.9 (Menard et al. 2003; Hamana et al. in prepa- 
ration; TJ03c). The iV-body simulations were carried out by the 
Virgo Consortium 8 (see also Yoshida, Sheth & Diaferio 2001), and 
were run using the particle-particle/particle-mesh (P 3 M) code with 
a force softening length of l so f t ~ 30/i~ kpc. The initial matter 
power spectrum was computed using CMBFAST (Seljak & Zal- 
darriaga 1996). For the analytic model, we approximate the initial 
condition to use the CDM transfer function given by Bardeen et 
al. (1986) with the shape parameter in Sugiyama (1995) for sim- 
plicity. The iV-body simulation employs 512 3 CDM particles in a 
cubic box of 479/i _1 Mpc on a side, and the particle mass of the 
simulation is m par t = 6.8 x 10 10 ft -1 Mq. 

The multiple-lens plane algorithm to simulate the lensing 
maps from the iV-body simulations is detailed in Jain et al. (2000) 
and Hamana & Mellier (2001). We will use the output data for 
source redshifts of z s = 1 and 3 for the following analysis. The 
simulated map is given on 1024 2 grids of a size of # gr id = 0'.2; the 
area is Q B = 11.7 degree 2 . The angular resolution that is unlikely 
affected by the discreteness of the iV-body simulation is around 1' 
(Menard et al. 2003; TJ03c). 

8 see 

http : / / star-www . dur .ac.uk/~frazerp/virgo/virgo. html 
for the details 
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Figure 2. The probability density function (PDF) of the magnification mea- 
sured from the ray-tracing simulations for z s = 1 and 3. To clarify the 
angular resolution, we show the two results of using different smoothing 
scales, which were used to avoid the artificial discrete effect of the N- 
body simulations. The comparison manifests that high magnification events 
are sensitive to smaller structures that are relevant for the smoothing. The 
dashed line shows an asymptotic behavior of the PDF for high magnification 
events, PDF oc 2 , as theoretically expected (see text in more detail). 

5.2 Probability distribution function of the magnification: 
the angular resolution of the simulations 

Universal properties of the critical curve, as demonstrated in §4, 
lead to an asymptotic behavior of the probability distribution func- 
tion (PDF) of large magnification events, irrespective of details of 
the lensing mass distribution: PDF(/i) oc fi~ 2 (Peacock 1982; 
Vierti & Ostriker 1983; Blandford & Narayan 1986; Schneider 
1987; Hamana, Martel & Futamase 2000, and also see Chapters 6, 
1 1 and 12 of Schneider et al. 1992). Note that the PDF is defined in 
the image plane, while the PDF becomes fj,~ 3 if one defines it from 
the sample of sources in the source plane. We use this property to 
investigate the angular resolution of the ray-tracing simulation. 

Figures 2 shows the magnification PDF measured in the sim- 
ulations for source redshifts of z s = 1 and 3, respectively. To com- 
pute the PDF, we accumulate the counts in a given bin of /i from 
36 realizations and then normalize the PDF amplitude to satisfy 
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Jdfj, PDF(/i) = 1 over the range of ft measured. The PDF has a 
skewed distribution: most events lie in demagnification of fi < 1 
and rare events have high magnifications with \i 1 having a long 
tail. These reflect an asymmetric mass distribution in the large-scale 
structure as expected from the CDM scenario - the underdense re- 
gion can be seen preferentially in the void region with a typical size 
~ lOMpc, while the highly non-linear structures appear in dark 
matter halos on scales ,$ Mpc. As can be seen, the simulation of 
z 3 — 3 displays more pronounced evidence of asymptotic depen- 
dence PDF(/i) oc /i~ 2 for high magnifications (p, £ 10) than the 
result for z 3 = 1. 

To more explicitly clarify the resolution issue of the simula- 
tions, we show the two results of the different smoothing scales 
that were used in making the projected density field to suppress the 
discreteness effect of the iV-body simulations (see Menard et al. 
2003 and Hamana et al. 2003 for more details). They are named 
as 'small-scale smoothing' (solid curve) and 'large-scale smooth- 
ing' (dotted curve), respectively. The former is expected to have the 
angular resolution around 1' as discussed in Menard et al. (2003) 
and in TJ03c, while the latter employs a smoothing scale two times 
larger than the former. The effect of the large-scale smoothing is 
that it more smoothes out smaller scale structures of the mass dis- 
tribution that are resolved by the small-scale smoothing simulation. 
The comparison manifests that occurrence of high magnification 
events (n > 1) is very sensitive to the small-scale structures. For 
this reason, we will employ the small-smoothing simulation in the 
following, since our interest is to clarify the non-linear magnifica- 
tion effect on the magnification statistics. This result also implies 
that simulations with higher resolution could further alter the PDF 
shape especially at high magnification tail. 

5.3 The two-point correlation function of lensing 
magnification 

We now turn to investigation of the magnification 2PCF, £ M = 
(<5/i • Sfi), as it is possible to observe from size fluctuations on dis- 
tant galaxy images (see §2.2). Figure 3 shows the comparison of 
the halo model prediction (solid curve) with the measurement from 
simulations (triangle symbol) for source redshift z s — 1. Note that 
the error bar in each bin denotes the sample variance for a simu- 
lated area of 11.7 degree 2 , which is computed from 36 realizations, 
and the errors in neighboring bins are highly correlated. In this and 
following results, we mainly employ the maximum magnification 
cutoff /Umax = 8 in the halo model prediction as well as in the 
simulation result. If we ignore the shear contribution to the mag- 
nification (1), this cutoff corresponds to /i w 1 + 2k — 2.3 for 
the weak lensing approximation. The cutoff value is chosen so that 
strong lensing events are removed from the analysis, because such 
events likely have greater magnification /i ^ 10 (private communi- 
cation with M. Oguri). Figure 2 shows that this cutoff leads us to 
exclude the events in a high magnification tail of the PDF. 

Figure 3 shows that the halo model prediction well matches 
the simulation result. The 1-halo term provides dominant contri- 
bution to the total power on small scales ^ 3', while the 2-halo 
term eventually captures the larger scale signal (see Figure Al of 
TJ03c). It is worth noting that the shear field in /i (see equation (1)) 
contributes to the 2PCF amplitude by ~ 10% over the scales con- 
sidered. To make clear the importance of the non-linear magnifica- 
tion contribution (5/i ^ 1), the dashed curve and the square symbol 
are the halo model prediction and the simulation result for the weak 
lensing approximation £ M « 4{kk) = 4£ K . For this case, 4£ re can 
be also computed from the fitting formula of the non-linear mass 
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Figure 3. The two-point correlation function of the magnification field 
against the separation angle for the ACDM model and source redshift 
z s = 1. The solid curve shows the halo model prediction, while the trian- 
gle symbol is the simulation result. In most part of this paper, the maximum 
magnification cutoff £i ma x = 8 is employed (see text for the details). The 
error bar denotes the sample variance for a simulated area of 11.7 degree 2 . 
For comparison, the weak lensing predictions, leading to § M « A(kk), are 
shown by the halo model (dashed curve), the simulation (square symbol) 
and the fitting formula of Smith et al. (2003; dot-dashed curve). Note that 
the simulation result is slightly shifted in the horizontal direction for illus- 
trative purpose. The lower panel explicitly plots the contribution of non- 
linear magnifications, £/i/4£ K — 1, for the halo model prediction and the 
simulation result. 

power spectrum recently proposed by Smith et al. (2003; hereafter 
Smith03), which demonstrates another test of the accuracy of the 
halo model as well as of the simulation. 

The lower panel explicitly plots the relative difference, 
C^/4£ K — !• The simulation result is computed from the mean val- 
ues of £ M and 4£ K , and we do not plot the large error bar for il- 
lustrative purpose. The correction of the non-linear magnification 
amounts to 30% at 9 2', and the non-negligible contribution 
of ~ 10% still remains even on large scales £ 10'. Our model of 
the 2-halo term (18) correctly captures the non-linear effect seen 
in the simulation on the large scales. This large-scale correction 
is somewhat surprising, since it is naively expected that the weak 
lensing approximation is valid at these scales. Menard et al. (2003) 
also showed that the non-linear correction on the large scales can be 
fairly explained taking into account the higher-order terms 0(k 2 ) 
in the Taylor expansion of [i, though the method ceases to be ac- 
curate on small scales £ 5'. One advantage of the halo model is 
that it allows us to explicitly introduce the maximum magnifica- 
tion cutoff in the calculation, which allows a fair comparison with 
the simulation and probably with the actual observation. In other 
words, the results shown depend upon the cutoff value employed. 
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Figure 4. As in the previous figure, but for z s = 3. For higher source 
redshifts, the non-linear magnification effect becomes more significant as 
expected. 



If we use the cutoff values of /i max = 2 and 100, the deviation, 
£m/4£k - 1, becomes ~ 20% and ~ 40% at 9 = 1', respectively 
(also see Figure 5). 

Figure 4 shows the result for z s = 3, as in the previous fig- 
ure. It is clear that the non-linear magnification contribution leads 
to significant enhancement in the amplitude of the magnification 
correlation relative to the weak lensing prediction. The enhance- 
ment is ^ 40% at 9 ,$ 2'. The comparison with the previous fig- 
ure manifests that sources at higher redshifts are more affected by 
the non-linear magnification, as the sources have more chance to 
encounter intervening halos. It is worth noting that, if we do not 
apply any masking of high magnification events in the simulation, 
the statistical error in each bin becomes very large, which indicates 
the presence of events with /i ^> 1 in some realizations as shown 
in Figure 2. 



5.4 A quantitative test of the asymptotic behavior of the 
magnification statistics 

In what follows we quantitatively test the asymptotic dependence 
of the magnification statistics on large magnifications (/i ^> 1), as 
derived in §4. For this purpose, we consider three cases of p = 
0.5, 1 and 1.5 for the magnification 2PCF parameterized as £ = 
{8/j, p ■ S/j, f ), with source redshift z s — 3. Figure 5 shows how the 
2PCF amplitude depends on the magnification cutoff /i max used in 
the halo model predictions and the simulation result. The separation 
angle of 9 = 1' is considered and the curves are normalized by the 
predictions from the weak lensing approximation. The scale 9 = 1' 
is chosen based on the fact that the scale is in the non-linear regime 
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Figure 5. Dependences of the magnification 2PCF amplitude on the max- 
imum magnification cutoff used in the evaluations. For the magnification 
correlation parameterized as {5fi p ■ Sfi p ), the three panels show the results 
for p = 1, 0.5 and 1.5, respectively. The source redshift z s = 3 and the 
separation angle 9 = 1' are considered. All the curves are normalized by 
the predictions derived from the weak lensing approximation (fi to 1+2k). 
As discussed in §4, it is expected that the 2PCF amplitude diverges for 
p > 1 for the limiting case /J max — * oo, while it remains finite for p < 1 
(see equation (24)). This is verified by the halo model as well as by the 
simulation. 



and unlikely affected by the angular resolution of the simulation 
(TI03c). First, one can see that even the most conservative choice of 
(i max = 2 leads to decent difference between the correct treatment 
and the weak lensing approximation. The consequence derived in 
§4 is that the 2PCF amplitudes for p = 0.5, 1.0 and 1.5 have the 
dependences on ^ max given as £ oc (i^; 5 , ln(^ max ) and /j^ x 
for /x max 3> 1, respectively. It is obvious that this consequence is 
verified by the halo model prediction as well as by the simulation 
result for ^ max ~ 10. Although the halo model results for a = 
1.0 and 1.5 display a bend at ^ max ~ 30, we found that this is 
due to high magnifications between the radial and tangential critical 
curves in NFW halos. Most importantly, the 2PCF amplitude for 
p = 0.5 has a well converged value for fi m&K ^ 5: the amplitude 
changes by less than 10% over /i max = [10, 1000]. Therefore, the 
statistics with p < 1 have a practically great advantage because 
it is little affected by the uncertainty in specifying the maximum 
magnification cutoff in the analysis. 

Finally, we note that the results shown above are unchanged 
even if we consider the cross-correlation {8fi p ■ S 2D ), where 8 2D 
is the projected density fluctuation field (e.g., see equation (25)), 
because the asymptotic behavior is determined by the power of n 
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entering into the general correlation function (p p f), as derived in 
§4. 

5.5 Application to QSO-galaxy cross correlation 

In this subsection, we consider an application of the halo model 
to the QSO-galaxy cross correlation. The angular fluctuation field 
of galaxies on the sky is a projection of the 3D galaxy fluctuation 
field 8 3 along the line-of-sight, weighted with the redshift selection 
function f g (z) of the galaxy sample: 



8n 9 (6) 



-f 

Jo 



dzf g (z) S 9 (d A (z)0,z), 



where f g {z) is normalized as J^°dzf g (i 
paper, we employ 



(25) 



I. Throughout this 



f g (z)dz : 



(3z 2 



sgr(3//J) 



exp 



(0 



dz, 



(26) 



with P — 1.5 and zo = 0.3. This model leads to the mean red- 
shift as z mcan = Jdz zfg(z) = 0.45 and roughly reproduces the 
actual distribution of galaxies in the redshift galaxy catalog (e.g., 
Dodelson et al. 2002). 

However, the galaxy fluctuation field 8 s is not straightforward 
to model, since the galaxy formation is affected by complex astro- 
physical processes in addition to the gravitational effect. Recently, 
Jain et al. (2003) developed a sophisticated description of the mag- 
nification correlations based on the halo model as well as the semi- 
analytic galaxy formation model. In particular, it was shown that it 
is crucial to account for a realistic model to describe how galaxies 
populate their parent halo, the so-called halo occupation number, to 
make the accurate model predictions on arcminute angular scales. 
The halo occupation number strongly depends on types of galax- 
ies such red or blue galaxies. We here address how the non-linear 
magnification further modifies the model prediction. 

Before going to this study, we consider the cross-correlation 
between the magnification field and the dark matter distribution, 
which corresponds to an unrealistic case that the galaxy distribu- 
tion exactly traces the underlying mass distribution; 8 g = 8. This 
investigation is aimed at clarifying how the non-linear magnifica- 
tion effect remains after inclusion of a realistic model of the galaxy 
clustering, from the comparison of the results with and without the 
galaxy bias model. In addition, in this case we can compare the 
model prediction with the simulation result that is computed from 
the same A-body simulation we have used. Extending the method 
presented in §3 leads to the 1-halo term contribution to the cross- 
correlation: 

j2t 



XH , <TV dz 



d X 



dM n(M) 



M 
Po 



poo p2ir 

x / #/ d V ><l>[nl?'- 1 (\4> + 0\)-i\ E M (0), (27) 
J o Jo 

where Em (a;) is the normalized projected density of the NFW pro- 
file given by equation (26) in TJ03b. Similarly, one can derive the 
2-halo term of ^s, as done by equation (18). 

In Figure 6 we shows the results for source redshifts of z s — 1 
(left panel) and 3 (right panel), as in Figures 3 and 4. Note that both 
the results employ the same redshift selection function (26) to ob- 
tain the projected density field. We simply assumed a special case 
of s — 4/5 for the magnitude slope for the unlensed QSO number 
count. From the comparison with Figures 3 and 4, it is clear that 
the non-linear magnification contribution is weakened, due to the 



single power of fi entering into the two-point correlation compared 
to the magnification 2PCF. Nevertheless, it should be stressed that 
the non-linear correction has significant contributions of > 15% 
and 25% at 6 ,$ 1' for z s — 1 and 3, respectively. 

Next, we consider a model of the QSO-galaxy correlation that 
takes into account both the galaxy bias and the non-linear mag- 
nification effect. To do this, we use the halo occupation number 
(iVg(M)) to describe how many galaxies populate their parent halo 
of a given mass M, in an average sense (e.g., Seljak 2000; Guzik & 
Seljak 2002; Jain et al. 2003; TJ03b; Cooray & Sheth 2002). Sim- 
ply replacing M/po in equation (27) with (N g (M)) /n ga ,\ leads to 
the 1-halo term of the QSO-galaxy cross correlation: 



tiiv) = [ 

Jo 



d X-r-l^fs(z)-r dMn{My 



dxdVt dx 



"gal 



^[mm^W-I] £m(|</> + 0|), (28) 



where n ga i is the average number density at epoch z defined as 
Mgai = JdM n(M){N g (M)}. The cross-correlation thus depends 
on the first moment of {N g (M)} . Note that, on the other hand, 
the two-point correlation of galaxies depends on the second mo- 
ment, and it contains somehow uncertainty in modeling the sub- 
Poissonian process in the regime of (N g (M)) < 1. 

As stressed in Guzik & Seljak (2002) and Jain et al. (2003), it 
is probably accurate to assume that one of the (N s ) galaxies in a 
halo sits at the halo center and this has decent impact on the model 
predictions. On the other hand, we assume that the other ((N g ) — 1) 
galaxies follow the dark matter distribution within the halo. Fol- 
lowing the method of Jain et al. (2003), the part of the integrand 
function in the 1-halo term (28) can be expressed as 

poo r-'lit 

(N g (M)) / d4> *p<f> [mm s_1 (0) - i] £m(|4> + e\) 

Jo Jo 

= x- 2 [a4 5s-1 (0)-i] 

poo p'2ir 

+ / # / dk P 4>[i*M" 1 {4>)-i\ x M (\4> + 0\) 

Jo Jo 

x[{N g (M))-l], (29) 
for (iV g (M)) > land 

= [mm 3_1 W-1] (N g (M)), (30) 

for (N g (M)) < 1. Substituting these equations into equation 
(28) leads to the halo model prediction for the QSO-galaxy cross- 
correlation. 

To complete the model prediction, we need an adequately ac- 
curate model of the halo occupation number (N g (M)). We em- 
ployed the model in Jain et al. (2003), which was derived from the 
GIF A^-body simulations, coupled to a semi-analytic galaxy for- 
mation model (Kauffmann et al. 1999). The simulation result of 
(N g (M)) is well fitted by the functional form, 



<A 9 (M)> 



M 
Mo 



+ ,4cxp [-A,(log 10 (M)-M B ) 2 ] .(31) 



The parameter values are taken from Table 1 labeled as 'z = 0.06' 
in Jain et al. (2003), which reproduces the measurements for to- 
tal (blue plus red) galaxies at 2 = 0.06 in the GIF simula- 
tions. In the following, we employ a lower mass cutoff of M > 



This is also the case for galaxy-galaxy lensing as shown in Guzik & Sel- 
jak (2002) 
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Figure 6. The cross correlation function between the magnification field and the projected density fluctuation field for z s = 1 (left panel) and 3 (right panel), 
as in Figure 3. To get the projected density field, we assume the redshift selection function given by equation (26). 



10 11 h~ x M® and ignore the redshift evolution of (N S (M)) for 
simplicity. This model leads to the galaxy bias parameter 6 ga i = 
(l/fig a i) / dMn(M)b(M) (N g (M)) = 1.2 at z = in the large- 
scale limit, thus reflecting the fact that the modeled galaxies are 
biased objects relative to the dark matter distribution. 

Figure 7 shows the model predictions for the QSO-galaxy 
cross-correlation, as in the previous figure. For comparison, the 
dot-dashed curve is the result for the cross-correlation between the 
magnification and the projected dark matter distribution in the pre- 
vious figure. The comparison of the solid and dot-dashed curves 
manifests that the realistic model of galaxy bias largely modifies 
the cross-correlation, and the galaxy bias cannot be described by a 
simple linear bias on the small angular scales (Jain et al. 2003). It is 
also clear that the non-linear magnification correction is 10 — 25% 
on arcminute scales. Therefore, this result implies that an inclusion 
of the non-linear effect will be necessary to make an unbiased in- 
terpretation of the precise measurement expected from forthcoming 
massive surveys such as the CFHT Legacy Survey and the SDSS. 



5.6 Sensitivity of the magnification statistics to the halo 
profile properties 

As discussed above, one of the useful cosmological information 
extracted from the QSO-galaxy correlation measurement is infor- 
mation on the halo occupation number of galaxies, which in turn 
provides a clue to understanding of galaxy formation in connec- 
tion to the dark matter halo properties (see Jain et al. 2003 for the 
details). We here demonstrate another possibility of the magnifi- 
cation statistics (especially measured via galaxy size fluctuations) 



to address questions: what can we learn from the measurements? 
How is this method complementary to the established cosmic shear 
that probes the correlations of the convergence or shear fields (k or 
7)? To examine this, we focus on the non-linear relation between 
the magnification and the cosmic shear fields, as given by equation 
(1). The non-linear effect is more pronounced on smaller scales, as 
have so far been shown. Future massive surveys promise to mea- 
sure the magnification statistics as well as the cosmic shear even on 
sub-arcminute scales (Jain 2002; TJ03c; Jain et al. 2003). Within 
a picture of the halo model, the sub-arcminute correlation function 
is quite sensitive to the halo profile properties (TJ03b,c) and the 
measurement can be potentially used to constrain the properties, 
if the systematics is well under control. Hence, we here investi- 
gate the dependence of the magnification 2PCF on the halo profile 
parameters, the halo concentration and the inner slope of the gen- 
eralized NFW profile. These parameters are still uncertain obser- 
vationally and theoretically and have information on the dark mat- 
ter nature as well as properties of highly non-linear gravitational 
clustering on IMpc. Following TJ03c, we consider the param- 
eterization given as c(M,z) = r v i r /r s — co(Af/M*) -013 and 
p(r) oc r~ a (l + r/r s )~ 3+a , respectively. Our fiducial model so 
far used is given by (co, a) = (9.0, 1.0). For cases a — 0, 1 and 
2 we can derive analytic expressions for the convergence and shear 
profiles from which we can also compute the magnification profile 
(the expressions of the convergence fields are given in Appendix B 
in TJ03c). Note that in what follows we employ the virial bound- 
ary condition. The relevant angular scales are below the angular 
resolution of A-body simulations we have used. 

The left panel of Figure 8 shows the halo model prediction 
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Figure 7. The halo model predictions for the angular cross-correlation between the magnification and the foreground galaxy distribution, where the galaxy 
clustering is modeled by the halo occupation number (31). For comparison, the dot-dashed curve is the result for the cross-correlation between the magnification 
and the projected dark matter distribution as shown in the previous figure. The non-linear magnification correction remains to enhance the cross-correlation 
amplitude relative to the weak lensing approximation, even if the realistic model of the galaxy clustering is included. 



for the magnification 2PCF with varying the halo concentration, 
while the right panel shows the results with varying a. Here we 
consider (<5/i ' 5 • Sfi 0,5 ), because it is less sensitive to high magni- 
fication events (p 3> 1) and therefore observationally more robust 
(see Figure 5). In the weak lensing limit, the correlation can be 
approximated by the convergence 2PCF (dashed curves), which is 
measured by the cosmic shear measurement because SfjP 5 ~ k. 
Therefore, the difference between the solid and dashed curves re- 
flects contribution from the non-linear magnifications S/j, ^ 1. Ha- 
los with masses M > 1O 13 M0 provide dominant contribution of 
^ 80% to the total power over a range of non-linear scales 0'.l-3' 
(e.g., see Figure 14 in TJ03c). One can see that the magnification 
2PCF has stronger sensitivity to the halo concentration and depends 
on the inner slope in a different way from the convergence 2PCF. 
In TJOc, it was pointed out that the cosmic shear measurement in- 
troduces a degeneracy in determining these halo profile parameters, 
even provided the accurate measurement (see Figures 16 and 17 in 
TJ03c). The results in Figure 8 thus indicate that a joint measure- 
ment of the magnification statistics and the cosmic shear can be 
used to improve the parameter determinations. Finally, one caution 
we make is that the magnification 2PCF for the profile with a — 2 
is more amplified by an increase of the maximum magnification 
cutoff /i m ax than the other a's and thus is sensitive to the selection 
effect. 



6 DISCUSSIONS 

In this paper, we have used the real-space halo approach to compute 
the magnification correlation function without employing the weak 
lensing approximation fi ~ 1 + 2k + 0(k 2 ). It has been shown 
that the correction due to the non-linear magnification (Sfi ^ 1) 
leads to significant enhancement in the correlation amplitude rela- 
tive to the weak lensing approximation (see Figures 3-7). The cor- 
rection is more important as one considers the correlation function 
for sources at higher redshifts and on smaller angular scales, where 
the weak lensing approximation ceases to be accurate. Thus, ac- 
counting for the non-linear contribution in the theoretical model is 
needed to extract unbiased, cosmological information from the pre- 
cise measurement expected from forthcoming and future surveys. 
The encouraging result shown is the halo model prediction remark- 
ably well reproduces the simulation result over the angular scales 
we consider. 

We also developed the model to predict the QSO-galaxy cross- 
correlation by incorporating the realistic model of the halo occupa- 
tion number of galaxies into the halo model (see §5.5). The pri- 
mary cosmological information provided from the measurement 
is constraints on the halo occupation number, as shown in Jain et 
al. (2003; also see Guzik & Seljak 2002). In particular, the QSO- 
galaxy correlation can be used to directly constrain the first moment 
of the halo occupation number, compared to the two-point correla- 
tion of galaxies that probes the second moment. Exploring the halo 
occupation is compelling in that it provides useful information on 
the galaxy formation and the merging history in connection with the 
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Figure 8. The dependences of the magnification 2PCF, {<5^ 0,5 ■ S[i 0,5 ), on the halo profile parameters. The left panel shows the halo model predictions for the 
halo concentration parameter co = 3, 9 and 15 (from bottom to top), while the right panel shows the results with a = 0, 1 and 2 for the inner slope parameter 
of the generalized NFW profile (see text for more details). For comparison, the dashed curves are the corresponding results for the convergence 2PCF. Note 
that the magnification 2PCF considered becomes equivalent to the convergence 2PCF in the weak lensing limit as S/J, ~ k. The lower panel shows the 
difference relative to the result with our fiducial model of a = 1 and co = 9. 



dark matter halo properties. We showed that the non-linear mag- 
nification amplifies the cross-correlation amplitude by 10 — 25% 
on arcminute scales. The method of this paper therefore provides 
the accurate model prediction that accounts for both the non-linear 
magnification correction and the realistic galaxy bias. 

We found that the magnification statistics can be used to ex- 
tract cosmological information complementary to that provided 
from the cosmic shear measurement. We have demonstrated that 
the joint measurement on angular scales 3' could be used to 
precisely constrain the halo profile properties (see Figure 8). This 
possibility would open a new direction in using the magnification 
statistics as a cosmological probe beyond determination of funda- 
mental cosmological parameters (Bartelmann 1995; Bartelmann & 
Schneider 2001; Menard & Bartelmann 2002; Menard et al. 2003). 
Exploring the halo profile properties with gravitational lensing will 
be a direct test of the CDM scenario in the highly non-linear regime 
^ Mpc, since alternative scenarios have been proposed in order 
to reconcile the possible conflicts between the CDM predictions 
and the observations on the small scales (e.g., Spergel & Steinhardt 
2000). 

In most results shown, we employed the maximum magnifi- 
cation cutoff (i m „ = 8 for the halo model predictions as well as 
for the simulation results, because the choice likely removes strong 
lensing events (pi £ 10) from the analysis. Even if we employ the 
smaller value, the qualitative conclusions derived are not largely 
changed, as can be seen from Figure 5. Observationally, strong 
lensing event is easily removed from the sample of the magnifi- 



cation statistics, if it accompanies multiple images or largely de- 
formed images. However, without the clear signature, it is relatively 
difficult to make a clear discrimination of the strong lensing, since 
the magnification is not a direct observable. One advantage of the 
halo model developed in this paper is that it allows a fair compari- 
son with the measurement by employing the selection criteria in the 
measurement for the model prediction. Based on these considera- 
tions, we derived useful, general dependences of the magnification 
correlation amplitude on large magnifications (/i 3> 1), from the 
universal lensing properties in the vicinity of critical curves (see 
§4). The intriguing consequence is that, for a correlation function 
parameterized as (fJ, p f), the amplitude converges to be finite for 
p < 1 and otherwise diverges p > 1 as the maximum magnifica- 
tion cutoff ptmax — > oo, independent of details of the lensing mass 
distribution. This was quantitatively verified by the halo model pre- 
diction as well as by the simulations (see Figure 5). This result 
therefore implies that the magnification statistics with p < 1 are 
practically advantageous in that it is insensitive to the selection ef- 
fect of the magnification cutoff /i ma x- This is the case for the two- 
point correlation of size (not area) fluctuations of distant galaxy 
images and for the QSO-galaxy cross correlation, if the unlensed 
number counts of QSOs with a limiting magnitude have a slope of 
s = d\aN{m)/dm < 4/5. 

One might imagine that the non-linear magnification contri- 
bution can be suppressed by clipping regions of cluster of galaxies 
from survey data in order to avoid the uncertainty in the model pre- 
diction and to apply the weak lensing approximation (e.g., see Bar- 
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ber & Taylor 2002). However, this likely adds an artificial selection 
effect in the analysis and causes a biased cosmological interpreta- 
tion. In addition, the lensing projection makes it relatively difficult 
to correctly identify the cluster region from the reconstructed con- 
vergence map, unless accurate photo-z information or follow-up 
observations are available (e.g., White et al. 2002). The approach 
of this paper allows us to treat data more objectively. 

There are some effects we have so far ignored in the halo 
model calculation. Most important one is the asphericity of halo 
profile in a statistical sense. The aspherical effect could lead to 
substantial enhancement of the magnification correlation ampli- 
tude, since it is known that an area enclosed by the critical curve 
is largely increased by ellipticity of the lensing mass distribution, 
thus increasing the cross section of high magnifications to the cor- 
relation evaluation. For example, the number of strong lensing arcs 
due to clusters of galaxies is amplified by an order of magnitudes if 
one considers an elliptical lens model instead of an axially symmet- 
ric profile (e.g., Meneghetti et al. 2003; Oguri et al. 2003). For the 
same reason, substructures within a halo could have strong impact 
on the magnification correlation, as they naturally emerges in the 
CDM simulations. Hence, it is of great interest to investigate the 
effect on the magnification statistics with higher resolution simula- 
tions. 
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